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— Double, Mesa Edited by Sandman on August 11, 1977 2:69 PM 

DIRECTORY 

DoubleDefs: FROM "doubledefs" . 

InlineDefs: FROM 'Mnl inedef s" , 

StringDefs: FROM "StringDef s"; 

DEFINITIONS FROM InlineDefs. DoubleDefs; 

Double: PROGRAM IMPORTS StringDefs EXPORTS DoubleDefs = 
PUBLIC BEGIN 

DDivide: PROCEDURE [num, den: LongCARDINAL] 
RETURNS [quotient, remainder: LongCARDINAL] » 
BEGIN 

qq, count: CARDINAL; 
ITemp: LongCARDINAL; 
IF den.highbits = THEN 
BEGIN 
[quotient. highbits, qq] ♦- 

LongDivMod[[lowbits :num.highbits,highbits:0],den. lowbits] ; 
[quotient . lowbits , remainder . lowbits] <- 

Lon gD i vMod[[ 1 owb i ts : num. lowbits,highbits:qq], den. lowbits]; 
remainder . highbits «- 0; 
END 
ELSE 
BEGIN 
count ♦- 0; 

quotient, highbits <- 0; 
ITemp ^ den; 

WHILE ITemp. highbits # DO -- normalize 
ITemp. lowbits ♦• 

BITSHIFT[lTemp. lowbits, -1] + BITSHIFT[lTemp. highbits, 15]; 
ITemp. highbits *- BITSHIFT[lTemp. highbits. -1] ; 
count ♦- count + 1; 
ENDLOOP; 
qq ♦- LongOiv['"um, ITemp . lowbits] ; -- trial quotient 
qq <- BITSHIFT[qq, -count]; 

ITemp ♦- LongMult[den. lowbits ,qq]; -- multiply by trial quotient 
ITemp. highbits ♦- ITemp . highbits + den . highbi ts*qq; 
UNTIL DCompare[lTemp, num] ^ greater DO 

-- decrease quotient until product is small enough 
ITemp ♦- DSub[lTemp. den]; 
qq ♦- qq - 1; 
ENDLOOP: 
quotient. lowbits ^ qq; 
remainder ♦- DSub[num, ITemp] ; 
END; 
RETURN 
END; 

DMultiply: PROCEDURE [a,b: LongCARDINAL] 
RETURNS [product: LongCARDINAL] = 
BEGIN 

product ♦- LongMult[a. lowbits, b. lowbits]; 
product . highbits ♦- 

product. highbits + a. lowbits*b .highb i ts + a. highb i ts*b . lowbits ; 
RETURN 
END; 

DAdd: PROCEDURE [a.b: LongCARDINAL] RETURNS [LongCARDINAL] = 
BEGIN 

t: CARDINAL = a. lowbits; 
a. lowbits ♦- a. lowbits + b. lowbits; 
a. highbits <- a. highbits + b. highbits; 
IF a. lowbits < L THEN a. highbits «- a. h i ghb i ts + 1 ; 
RETURN[a] 
END; 

DSub: PROCrOURC [a,b: I ongCARDTNAL] RETURNS [LongCARDINAL] = 
BEGIN 

t: CARDINAL = a. lowbits; 
a. lowbits «- a. lowbits - b. lowbits: 
a. highbits ♦- a. highbits - b. highbits; 
ir a. lowbits > L THEN a. highbits <- a. h i ghb i ts-1 ; 
RrrURN[a] 
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END; 

DNeg: PROCEDURE [a: LongCARDINAL] RETURNS [LongCAROINAL] » 
BEGIN 

IF (a.lowbits ^ -a.lowbits) = THEN a.highbits ^ -a.highbits 
ELSE a.highbits ♦- BITNOT[a.highbits]; 
RETURN[a]; 
END; 

Dine: PROCEDURE [a: LongCARDINAL] RETURNS [LongCARDINAL] » 
BEGIN 
IF (a.lowbits *- a.lowbits + 1) = THEN 

a.highbits ♦- a.highbits + 1; 
RETURN[a] 
END; 

DCompare: PROCEDURE [a.b: LongCARDINAL] RETURNS [Comparison] « 
BEGIN 

IF a = b THEN RETURN[equal] ; 
RETURN[SELECT a.highbits FROM 
< b.highbits => less, 
> b.highbits -> greater, 
ENDCASE => 

IF a.lowbits < b.lowbits THEN less ELSE greater] 
END; 

AppendDouble: PROCEDURE [s: STRING, a: LongCARDINAL] = 
BEGIN OPEN StringDefs; 

longZero: LongCARDINAL = [highbits:0, lowbits:0]; 
longlO: LongCARDINAL = [highbits:0. lowbits:10]: 
xn: PROCEDURE = 
BEGIN 

charZero: CARDINAL « LOOPHOLE[ '0] ; 
r: LongCARDINAL; 
IF a # longZero THEN 
BEGIN 

[a, r] ♦■ DDivide[a, longlO]; 
xn[]: 

AppendChar[s, LOOPHOLE[r. lowbits+charZero, CHARACTER]]; 
END; 
END; 
IF a = longZero THEN AppendChar[s , '0] ELSE xn[]; 
RETURN 

END; — AppendDouble 
END.. . 



